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ABSTRACT 

We present observations of the singly and doubly deuterated species of thio- 
formaldehyde, HDCS and D2CS, towards the dark cloud Barnard 1. This is the 
first detection of D 2 CS in Space and in dense and cold prestellar regions. Col- 
umn densities obtained using rotational diagrams and a Large Velocity Gradient 
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model show an extremely high D-enhancement in thioformaldehyde in Barnard 1. 
Although the column density of H 2 CS is smaller than that of H 2 CO, both species 
show similar D-enhancements in their singly and doubly deuterated species. A 
chemical model -including multiply deuterated species- has been used in order to 
interpret the observations. 

Predicted rotational frequencies from laboratory data for HDCS and D 2 CS 
are significantly in error when compared to the observed frequencies in Space. 
Consequently, we have derived new rotational constants for these two species and 
for H 2 CS and H 2 C 34 S using the observed frequencies in Barnard 1. The new ro- 
tational constants allow to predict the rotational transitions of these species with 
the accuracy needed for the narrow line emerging from dark clouds. Rotational 
constants for HDCS and D 2 CS have been obtained from the observed transitions 
in the laboratory and in Space. 

Subject headings: ISM: abundances — ISM: individual objects: (Barnard 1) - 
ISM: molecules 



1. INTRODUCTION 

For a long time, singly deuterated interstellar molecules are known to be present in cold 
molecular clouds with an extremely high abundance enhancement due to chemical fractiona- 
tion (van Dishoeck et al. 1995; Minowa et al. 1997; Saito et al. 2000; Tine et al. 2000; Turner 
2001; Parise et al. 2002). The formation of deuterated molecules is preferred at low temper- 
atures (<80 K) and leads to a high degree of fractionation in cold and dense dark clouds 
without embedded luminous sources. These prestellar cores offer the opportunity to study 
the cold and dense evolutionary stage in which the deuteration of molecules supposedly oc- 
curred. Although this is a relatively well understood chemical process, the observation in the 
last few years of multiply deuterated molecules has renewed interest in deuterium fraction- 
ation in dense clouds : ND 2 H (Roueff et al. 2000); ND 2 H and D 2 CO (Loinard et al. 2001); 
ND 3 (Lis et al. 2002); D 2 S (Vastel et al. 2003); D 2 CO (Ceccarelli et al. 2002; Bacmann et 
al. 2003); CHD 2 OH and CD 3 OH (Parise et al. 2002, 2004). Among these deuterium isotopic 
species, there are, so far, only three sulfur bearing deuterated molecules, HDS and D 2 S, 
detected towards several protostars and dense cores (Vastel et al. 2003; van Dishoeck et al. 
1995); and HDCS, which has been detected toward TMC-1 (Minowa et al. 1997). 

Our new data and, in particular the high enhancement in the fractionation of multiply 
deuterated species, are requiring a review of chemical models. The presence of deuterated 
molecular species in star forming regions like Orion-IRc2 and IRAS 16293-2422 (Turner 
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1990; van Dishoeck et al. 1995; Ceccarelli et al. 1998; Loinard et al. 2000) has been taken as 
signs for an active grain surface chemistry followed by desorption when the physical processes 
associated to recent star formation activity increases the temperature of the gas and of dust 
grains. However, this is not the case for dark clouds, where the kinetic temperature is too 
low for efficient evaporation of ice grain mantles. Roberts & Millar (2000) found that gas 
phase chemistry can account for the large abundances of both singly and doubly deuterated 
molecules when the gas is very depleted of molecules and atoms. The reason is that both 
CO and O efficiently destroy Hg . Hence when CO and O are sufficiently depleted, Hij" is 
converted efficiently into H 2 D + in the reaction with HD. Subsequent reactions of H 2 D + with 
HD leads to D 2 H + and finally to the formation of Dg . These ions propagate the deuterium 
to other molecules (Lis et al. 2002; Roberts et al. 2003). In their model, Roberts et al. (2003), 
conclude that in regions where CO is higly depleted, the abundance of D 2 H + will be simillar 
than that of H 2 D + and that could be the most abundant form of deuterated Hj~. The 
first prediction has been confirmed recently with the detection of D 2 H + toward the prestellar 
core 16293E (Vastel et al. 2004). For an early evolution time (~2-4xl0 4 yr) the double and 
triple deuterated species could have large enhancements (Roberts et al. 2003). In the case 
of Barnard 1, Lis et al. (2002), based on their C 18 and C 17 (2-1) observations found that 
C 18 O/H 2 =5.4xl0~ 8 , which is ~3 times lower than the canonical value, C 18 0/H 2 =1.7x 10~ 7 
(Frerking et al. 1982; Cernicharo & Guelin 1987), indicating that an appreciable amount of 
CO is depleted onto grains. Modelling with low temperature conditions and high depletion 
factors allows a satisfactory comparison with the observations of doubly and triply deuterated 
ammonia -see e.g. NHD 2 towards L134N (Roueff et al. 2000)- and also helps to explain the 
high abundances observed in star forming regions in terms of grain chemistry models. The 
recent detection of doubly deuterated H 2 S provides a stimulating impetus to consider other 
sulfur molecules and extend the previously elaborated chemical schemes. 

During a 3mm line survey on dark clouds (Barnard 1, TMC-1, L1544, L183 -Marcelino 
et al., in prep.), we found several transitions that initially were considered as unidentified. 
Taking into account the reduced number of possible candidates, the fact that the lines were 
strongest in Barnard 1 (see Figure 1), and the strong deuterium enhancement found in these 
sources by several authors (see, e.g., Roueff et al. 2000; Gerin et al. 2001, and references 
therein), we concluded that these lines were arising from D 2 CS and HDCS. While single 
deuterated thioformaldehyde had been detected at lower frequencies in TMC-1 (Minowa et 
al. 1997); this is now the first detection of D 2 CS in Space. Once the carriers were identified, 
and given the extremely large deuterium enhacement indicated by the observations, we 
concentrated our efforts and available observing time in the source showing the strongest 
lines, i.e., Barnard 1. In Section 2 we present the observing procedures and in Section 3 the 
main results. In this section we discuss the rotational constants of thioformaldehyde and its 
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isotopomers. We also analyze the physical and chemical conditions of Barnard 1 and derive 
the column densities for the different thioformaldehyde isotopic species. Section 4 is devoted 
to the discussion and to the chemical analysis of the results. 

2. OBSERVATIONS 

The observations were carried out using the IRAM 30-m radio telescope at Pico Veleta 
(near Granada, Spain) between 2003 October and 2004 February. Spectral lines arising from 
singly and doubly deuterated thioformaldehyde between 83 and 155 GHz were observed in 
October. In January we performed most of the observations for the main isotopic species. 
Up to four SIS receivers were used simultaneously. Each receiver was tuned to single side 
band with image rejections within 20-27 dB (3mm receivers), 12-16 dB (2mm receivers) and 
13 dB (1mm receivers). System temperatures were 100-225 K in October for the 3 and 2mm 
receivers with atmospheric opacities between 0.05-0.1. In 2004 January weather conditions 
were even better allowing excellent observations at 1mm. During this run we obtained 
temperatures below 100 K at 3mm and below 200 K at 1mm (atmospheric opacities at 225 
GHz were below 0.1). 

Calibration was performed using two absorbers at different temperatures and the ATM 
code developed by Cernicharo (1985) and recently updated by Pardo et al. (2001). An 
autocorrelator providing several backends with 40 kHz and 80 kHz of resolution was used. 
The velocity resolution was between 0.08-0.23 kms -1 depending on the frequency. Telescope 
pointing and focus were checked every 1 and 3 hours respectively on strong and nearby 
continuum sources. The largest pointing error between consecutive pointing scans never 
exceeded 3". 

Most of the observations were performed using wobbler and frequency switching tech- 
niques, using for the latter a frequency throw of 7.2 MHz -which eliminates the main ripple 
component coming from the telescope- and in some cases we used position switching. The 
throw for wobbler and position switching was 240" and 10' respectively, both in azimut. 
Taking into account the size of the source, these offsets are enough to provide a clean refer- 
ence. On-source integration times were typically ~l-2h, resulting in rms main beam antenna 
temperatures between 15-50 mK. Observed line parameters obtained from Gaussian fits to 
the observed lines are given in Table 1 and the line profiles are shown in Fig. 2-4. Line inten- 
sities for the observed transitions are given in main beam brightness temperature scale. The 
main beam telescope efficiencies are 0.78-0.75, 0.69 and 0.57 at 3, 2 and 1mm respectively. 

In order to compare the deuterium enhancement in thioformaldehyde with that of 
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formaldehyde we also observed several rotational transitions of H 2 CO, HDCO, D 2 CO and 
H^CO in 2004 February. The observed line parameters are given in Table 2 and the corre- 
sponding line profiles are shown in Figure 5. 

3. RESULTS 

Figure 2a shows the observed lines for the singly and doubly deuterated species of 
thioformaldehyde. The fact that the same transitions of both molecules appear at different 
velocities when assuming the predicted frequencies from available line frequency catalogs 
-see, e.g., JPL (Pickett et al. 1998); Koln (Miiller et al. 2001); Lovas (Lovas 1992) and 
Cernicharo (Cernicharo et al. 2000)- indicates a considerable error in the frequencies used. 
We have used the new frequencies for these species derived from our observations to improve 
the molecular constants and to obtain better values for line prediction at higher frequencies 
(see Section 3.1). Figure 2b shows that the agreement between HDCS and D 2 CS line profiles 
is excellent when the new predicted frequencies for both species are used. 

In Figure 2 we can see that lines arising from HDCS are moderately stronger than those 
from D 2 CS, which, although indicating a large column density for the singly deuterated 
species also implies a very high D-enhancement for D 2 CS (see below). Line widths are 
rather narrow as indicated in Table 1. The similarity between line profiles from both species 
indicates that they coexist spatially. 

We observe a double peak profile for H 2 CO lines probably due to selfabsorption by 
foreground gas (the high velocity gas associated with IRAS 03301+3057, a young stellar 
object which is close to the strong ammonia peak, could also, perhaps, affect the line profile of 
formaldehyde; we have checked that the absorption feature is not coming from emission in the 
off position by checking the spectra in frequency switching mode). Lines from the deuterated 
species of H 2 CO and from H 2 13 CO are wider than the lines arising from thioformaldehyde and 
its isotopomers, but less than those of H 2 CO and they do not show the line wings observed 
in formaldehyde. Lines arising from H 2 CO are stronger than those from H 2 CS and lines 
from HDCO and D 2 CO are about two times stronger than the lines arising from deuterated 
thioformaldehyde, indicating a higher column density for the different H 2 CO isotopic species. 

3.1. Molecular Line Frequencies 

In order to have a coherent picture of the line profile from the different transitions 
of thioformaldehyde and its isotopes we performed a determination of the velocity of the 
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source as a mandatory step in deriving line frequencies. We used 19 molecular lines that 
were observed in Barnard 1 during the present work, and for which laboratory frequencies or 
confident frequency predictions are available (JPL, Koln, Lovas, and/or Cernicharo catalogs; 
see above). The lines showing some evidences of opacity effects were discarded. Finally, only 
lines from CCH (some hyperfine components), methanol, deuterated formaldehyde, CCS, 
and one of the hyperfine components of the N = 1 — transition of CN were used. The 
averaged velocity obtained from these lines is 6.70±0.05 km s _1 , where the error translates 
the uncertainty from laboratory measurements and also the presence of a second velocity 
component in the cloud at 7.5 km s _1 that reduces the accuracy of the Gaussian fits. The 
estimated LSR velocity of 6.70 km s _1 was assumed for all lines in Table 1 when deriving 
observed frequencies. 

We have observed four HDCS transitions, based on the frequencies calculated from the 
rotational constants of Minowa et al. (1997) and seven from D 2 CS based on the frequencies 
derived from the rotational constants of Cox et al. (1982). The rotational constants for the 
latter species have been derived from high J transitions. However, the largest frequency 
covered by the laboratory data is 35 GHz, hence, one could expect significant uncertainties 
for rotational transitions at large frequencies. A comparison of the line profiles of HDCS 
and D 2 CS (Fig. 2) suggests frequency errors in the latter species as large as 0.1-0.2 MHz. 
Even for HDCS for which the rotational constants have been derived from a significant set of 
data (Minowa et al. 1997), observed frequencies in dark clouds are often different from those 
obtained in the laboratory. As an example, the 3 3-2 02 line of HDCS has been observed by 
Minowa et al. (1997) with a laboratory frequency of 92981.601±0.018 MHz, while the same 
authors reported a frequency for this line from observations of TMC-1 of 92981.680 MHz. We 
have observed the same line in Barnard 1 with much better signal to noise ratio and we obtain 
(see Table 1) a frequency of 92981. 658±0. 030 MHz. Consequently, the use of the laboratory 
frequency will produce a 0.17 km s _1 velocity shift in the line -see Figure 2. It seems evident 
that the calculated frequencies in the available line catalogs have significant errors for some 
molecules. This is a critical point, because astronomers often compare line profiles from 
different isotopomers of a given molecular species to trace cloud collapse or line opacities. 
Uncertainties in computed frequencies larger than 20 kHz will produce significant velocity 
shifts when comparing different lines from the same species or from different molecules. 

The situation is even worse for rare isotopic substitutions like deuterium, 34 S or 13 C. The 
rare isotopomer H 2 C 34 S is a much better tracer of the column density of thioformaldehyde 
than the lines from the main isotope. We observed several lines of this species (see Table 1) 
using frequencies from the JPL line catalog. The number of measured laboratory frequencies 
for this species is small and, hence, computed frequencies have considerable uncertainties. 
For example the 3 3-2 2 line predicted at 101283.899 MHz was observed at 101284.331 MHz, 



-7- 



the 3i2-2n line predicted at 102810.609 MHz was observed at 102807.426 MHz, the 4 14 -3i 3 
predicted at 133031.473 MHz was observed at 133027.035 MHz, and the 4 13 -3i 2 line predicted 
at 137077.543 MHz was observed at 137071.099 MHz. That the observed frequencies were 
corresponding to the searched lines of H 2 C 34 S was confirmed by a simultaneous fit (see 
below) to the laboratory and space frequencies of a few rotational constants (several of them 
fixed to the values of the main isotope -note that H 2 CS and H 2 C 34 S have the same inertia 
moment around the A axis of the molecule). The large differences that are found between 
space spectroscopy and predicted frequencies could prompt for negative results and/or to 
reach a confuse situation as we often look the on-line frequency catalogs as a final product 
ready to be used. Special care has to be taken when using frequencies predicted from a fit 
to a reduced number of observed lines, or for predictions in frequency domains outside those 
observed in the laboratories. 

Using the astronomically measured frequencies given in Table 1 and those observed in 
the laboratories we performed fits to the rotational constants of H 2 CS, H 2 C 34 S, D 2 CS and 
HDCS. The results are given in Tables 3 to 6. These tables also give the fits obtained 
from the laboratory data alone. Although HDCS is a nearly prolate symmetric top and a 
S-reduced Hamiltonian should be preferred (Gordy & Cook 1984; Watson 1977; Minowa 
et al. 1997), we adopted the asymmetric (A) reduced form of the Watson Hamiltonian in 
a F representation (Watson 1967, 1977). Laboratory frequencies for the different isotopic 
species have been taken from Clouthier & Ramsay (1983), Dangoise et al. (1978), Chu et 
al. (1973), Cox et al. (1982) and McNaughton & Bruget (1993). Although it seems that 
the rotational constants in Tables 3-6 show little variation between those obtained from 
laboratory data and those derived from laboratory plus space observations, these variations 
account for changes in the predicted frequencies of 0.1-0.2 MHz between both sets of data. 
The species the less affected by the introduction of the space frequencies is HDCS and the 
most H 2 C 34 S. 

One could think that the situation should be better for the main isotope of thioformalde- 
hyde for which 60 of its rotational lines have been measured in the laboratory. However, 
Figure 3 shows that the laboratory frequencies and the calculated ones result in different ve- 
locities for the lines. Although, the uncertainties in the laboratory measurements are large, 
we could expect that the global fit to all laboratory data will produce rotational constants 
accurate enough to allow a good estimate of the spectrum of H 2 CS. As can be seen in the 
central panels of Figure 3 this is really the case, but these predictions are not accurate enough 
for the narrow lines emerging from dark clouds. Hence, the frequencies of the eight observed 
lines of H 2 CS were obtained following a similar procedure to that adopted for its rare iso- 
topes and a new set of rotational constants was obtained. The results are given in Table 3. 
With the new set of rotational constants for the four isotopomers of thioformaldehyde, we 
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have computed their frequencies and those corresponding to the observed lines are given in 
Table 1. 

The observed lines of H 2 CO show an asymmetric profile with a red line wing which could 
arise from molecular outflows (see below). For H 2 CO, HDCO and also HDCS, the observed 
lines present a second peak at ~7.5 km s _1 . This velocity component is also observed in 
HC 17 0+, CCH, c-C 3 H and HC 3 N. Bachiller & Cernicharo (1984) already found, based on 
the asymmetric and wide 13 CO line profiles, two velocity components in the cloud (at 6.5 
and 7.2 km s _1 ). The lines of D2CO peak at 6.7 km s -1 (see Figure 5), while those of HDCO 
and Hj, 3 CO have their emission maxima at slightly lower velocities. If we assume for these 
lines a cloud velocity of 6.7 kms -1 , then the frequencies for the lines of these two isotopes 
of H 2 CO will be 0.1 MHz above of the predicted ones, which are much more accurate (large 
number of measured lines in the laboratory) than those of the isotopes of H 2 CS. Hence the 
observed velocities for HDCO and Hj, 3 CO could correspond to optical depth effects or to 
chemical cloud structure. In fact, the high intensity of the HDCO lines and the asymmetric 
profiles in H2CO indicate high optical depths in the Barnard 1 core, or high foreground 
column densities, in the rotational lines of this species, even in those of HDCO. 



3.2. Column Densities and Physical Conditions 

3.2.1. Previous Studies of Barnard 1 

Barnard 1 is one of the highest column density molecular sources in the Perseus complex 
with recent low-mass star formation activity, as indicated by several young stellar objects - 
three of them without optical counterpart- detected by IRAS. Bachiller & Cernicharo (1984) 
made the first study of this region using visual extinction maps derived from star counts (see 
Cernicharo et al. 1984), and emission maps of CO and its isotopomers, HCO + , H 13 CO + 
and NH 3 . From star counts, they found Barnard 1 to be composed of two condensations of 
A„>5 mag. They modeled the cloud in four regions -core, internal and external envelope 
and halo- of different size and physical conditions depending on the line intensities of the 
observed molecules (see also Cernicharo & Guelin 1987). The innermost layer, i.e. the 
core, was studied from ammonia data. With a visual extinction ~10 and a size of 0.5 pc 
(at 200 pc of distance), the average density in the core ranges from 10 3 to 3xl0 4 cm -3 and 
the kinetic temperature is 12±2 K (in concordance with CO data). H 13 CO + and HCO + 
emission is also detected in this region with column densities N(H 13 CO + )~2xl0 12 cm~ 2 and 
N(HCO + )~2xl0 14 cm -2 , where the last one extends to the internal envelope. The rest of 
the cloud can be characterised by different isotopomers of CO. The density of the internal 
envelope was obtained from C 18 emission and results to be 2xl0 3 cm -3 . This region is 
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very extended (~1 pc) and contains the young stellar objects and the NH 3 globule. The 
visual extinction here is ~3-4 mag. The external envelope, with a density below 10 3 cm -3 
and A„~2-3 mag, was studied from 13 CO emission. The halo (A„~l-2 mag) is the most 
external zone with a large spatial size and traced by moderate/weak CO emission -see also 
Cernicharo & Guelin (1987). They also made a dynamical study of the cloud and concluded 
that it is being supported by magnetic fields and/or turbulence. 

Bachiller et al. (1990) mapped Barnard 1 in several molecules with higher angular 
resolution, resolving some spatial structures in the main dense core from observations of 
NH3. In this region they found density values of 7xl0 4 cm -3 in the southern clump, and 
>3.7xl0 4 cm -3 in the northern one. They studied the external and the internal envelopes 
from 13 CO, C 18 and CS emission. CS maps show the dense region to have a size of ~ 2 
pc x 5 pc (at 350 pc of distance). The strongest emission of this molecule is restricted to 
the internal envelope, but extends to the external one with lower intensities, probably due 
to the existence of several dense clumps. They obtained n(H 2 )>2xl0 3 and 10 3 cm -3 for the 
internal and the external envelopes, respectively. The kinetic temperature is 11-14 K in the 
densest regions and increases to Tk; n >20 K at the edge of the cloud, showing that the cloud 
is submitted to some heating from the exterior. 

Continuum observations indicate the presence of four clumps in Barnard 1 (see Matthews 
& Wilson 2002). Among them, Bl-b corresponds to the position we are observing and the 
source IRAS03301+3057 is close to Bl-a. Barnardl-b has two intense (sub) millimeter con- 
tinuum sources, probably Class sources associated to compact protostellar objects (Hirano 
et al. 1999). However, no outflows have been detected from these objects. Based on their 
dust continuum observations, Hirano et al. (1999) obtained a H 2 column density ~1.3xl0 23 
cm~ 2 in a 15" beam. Assuming a distance of 200 pc, this column density corresponds to a 
mean gas density of 2.9 xlO 6 cm~ 3 in the beam diameter, showing that high gas densities 
are present in the source. From CO J = 2 — 1 observations, Bachiller et al. (1990) found a 
high velocity flow emerging from IRAS 03301+3057, a young stellar object located near the 
strong ammonia peak (Bl-a). This outflow, which covers a region of about 40", was studied 
in detail by Hirano et al. (1997) from CO (J = 1 — 0) emission and by Yamamoto et al. 
(1992) from SiO J = 2 — 1 and J = 1 — lines. The interaction of the flow with the ambient 
gas is traced by the different Herbig-Haro objects present in this cloud (Yan et al. 1998). 

The wings in the H2CO lines shown in Figure 5 could be produced by the high velocity 
gas associated to IRAS 03301+3057 (located in Bl-a). Note, however, that our position 
(the Bl-b clump) is about 1' East from this source, and from the compact SiO peak which 
marks the place were the outflow interacts with the cloud material. Hence, the wings could 
be unrelated to the deuterium source. 



3.2.2. Molecular Background 



Formaldehyde and thioformaldehyde have both a nuclear spin statistics that separates 
the rotational levels into two different sets. Radiative and collisional transitions between 
ortho and para species are strictly forbidden, hence, these two species behave like two dif- 
ferent molecules. Rotational levels with k_ x and k +1 values even/odd or even/even belong 
to the para species while those with odd/even and odd/odd belong to the ortho species. For 
H 2 CO, H 2 CS, and their 13 C isotopic species the nuclear spin statistics leads to a 3:1 value 
for the ortho/para ratio. The singly deuterated species HDCO and HDCS do not have this 
separation on the rotational levels because the symmetry of the molecule has been broken 
by the deuteration. However, these species have b-type rotational transitions as deuteration 
introduces a weak permanent dipole moment along the B axis of the molecule (typically 

0. 2-0.3 D). Doubly deuterated species on the other hand also have a symmetry axis, and the 
nuclear spin statistics separate the rotational levels. However, as the nuclear spins are differ- 
ent, rotational levels of doubly deuterated isotopic species with k_ x and k +1 values odd/even 
or odd/odd belong to the para species, and those with even/odd even/even to the ortho one, 

1. e., just the opposite than for the main isotopic species. Moreover, the ortho/para ratio in 
the doubly deuterated species of formaldehyde and thioformaldehyde is 2. 



3. 2. 3. Thioformaldehyde 

In order to derive column densities for thioformaldehyde and their isotopes we have 
proceeded in two different ways. The first method we have used consists in the analy- 
sis of rotational diagrams for the observed isotopes of thioformaldehyde, see, e.g., Gold- 
smith & Langer (1999). From these diagrams one can obtain the molecular column den- 
sity and the excitation temperature of the molecular rotational levels (assuming all levels 
are populated by a single excitation temperature). From these diagrams we obtain com- 
mon rotational temperatures close to 10 K, T rot =9.5±1.0, 10.2±2.0, and 7.4±0.8 K for 
HDCS, D 2 CS and H 2 CS respectively, which suggests that the molecular hydrogen volume 
density is large enough to thermalize the rotational levels of these molecules (in the case 
of thermalization, T rot should be close to T kin ). The derived column densities for these 
species, assuming that the emitting source fills the beam of the 30-m radio telescope, are 
N(H 2 CS) = 1.2-2.0xl0 13 cm" 2 , N(HDCS)=4.5-7.0x 10 12 cm" 2 , N(D 2 CS) = 1.0-4.0x 10 12 cm" 2 . 
Hence, the abundance ratios are N(H 2 CS)/N(HDCS)~3±1, N(HDCS)/N(D 2 CS)~3±2, and 
N(H 2 CS)/N(D 2 CS)~9±3. Finally, for H 2 C 34 S we have assumed the rotational temperature 
obtained for H 2 CS to derive N(H 2 C 34 S)= 6.0±3.0xlO n cm" 2 . The column densities of H 2 CS 
and H 2 C 34 S are compatible with a galactic 32 S/ 34 S ratio of 24.4±5.0 (Chin et al. 1996). This 
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analysis does not provide any insight on the volume density and assumes that the lines 
are optically thin (which seems to be correct in view of the low column density derived for 
H 2 C 34 S). The volume density and other physical parameters could be derived through the 
observed line intensities and a Large Velocity Gradient code (see, e.g., Goldreich et al. 1974) 
if the collisional rates are known. We have adopted for H2CS the collisional rates of formalde- 
hyde computed by Green et al. (1978), multiplied by a factor of 2 to take into account the 
different effective surface and different energy of the levels of H 2 CS (a similar factor to that 
existing between the collisional rates of CO and of CS). Taking into account that we have 
observed a large number of lines of H 2 CS, we could check that the adopted collisional cross 
sections produce a good fit to the observed intensities using reasonable values for the gas 
density and temperature. Of course, the final result has to be in good agreement with the 
column densities and rotational temperatures derived from the rotational diagrams. Assum- 
ing a kinetic temperature of 10 K (see section 3.2.1) we derive rotational temperatures for 
the observed transitions of thioformaldehyde between 6.5 and 9.1 K for volume densities of 
2-4 xlO 5 cm -3 and a column density of ortho and para species of 1.2 and 0.5 xlO 13 cm~ 2 re- 
spectively. These values are in reasonable agreement with those derived from the rotational 
diagrams. 

For D 2 CS we have also used the same collision rates as for H 2 CS. Assuming also a kinetic 
temperature of 10 K, we derive n(H 2 )=3xl0 5 cm" 3 , T rot ~10 K, and N(D 2 CS)~3.0x 10 12 
cm -2 , with a derived value for the ortho/para ratio of ~ 1.7, very close to the expected value 
of 2. 

For HDCS the collisional rates could be rather different from those of the main isotope. 
We have tried to fit the data using the collisional rates for HNCO, but it was difficult to 
reach a good fit to the data with reasonable values for the physical parameters (which could 
be expected from this crude approximation to the collisional rates of HDCS). We adopt for 
HDCS the values obtained above from the rotational diagrams (note that the LVG values 
for the other species compare very well with those obtained from the rotational diagrams). 

3.2.4- Formaldehyde 

The column density of formaldehyde and its isotopes have been derived from an LVG 
analysis of the observed intensities, adopting the collisional rates of Green et al. (1978). The 
line profile of H 2 CO indicates opacities larger than 1. Hence, the column density for this 
species has been derived from the isotope H 2 3 CO assuming an 12 C/ 13 C isotopic ratio for dark 
clouds of 50 (Cernicharo & Guelin 1987; Langer et al. 1980). Assuming again a kinetic tem- 
perature of 10 K, the observed intensities give N(H 2 3 CO)=1.2±0.5xl0 12 cm" 2 , n(H 2 )=3xl0 5 
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cm -3 . There is no indication for an ortho/para ratio different than 3 from the observational 
data. The column density for the doubly deuterated species is N(D 2 CO)=3.0±1.5xl0 12 
cm -2 . For the single deuterated species we have assumed the same rotational temperatures as 
for the other isotopes of formaldehyde to obtain N(HDCO)=7±2xl0 12 cm~ 2 . Hence, the iso- 
topic ratios are N(H 2 CO)/N(HDCO)~7, N(H 2 CO)/N(D 2 CO)~13, and N(HDCO)/N(D 2 CO)~2.5. 
While the abundance ratios of singly over doubly deuterated species for formaldehyde and 
thio formaldehyde are similar (~3), the abundance ratio between the main isotope and the 
singly deuterated species is twice larger in thioformaldehyde compared to formaldehyde. 



4. ABUNDANCE RATIOS, DEUTERIUM FRACTIONATION AND 

CHEMICAL MODELLING 

Table 7 sumarizes several studies in dark clouds and star forming regions where deuter- 
ated molecules have been found. For Barnard 1, the new ratios obtained in this work for 
H 2 CS and H 2 CO are also shown. Comparing with other deuterated species detected towards 
this core, we can see that the D-enhancement of thioformaldehyde is among the highest 
values. In general the abundance ratios in Barnard 1 are only comparable with those found 
in the class binary protostar IRAS 16293-2422. That means that deuteration processes 
should be similar in both sources and different than in dark and quiescent clouds. 

We have extended our chemical network dedicated to deuterium containing molecules 
(Lis et al. 2002, Roueff et al. 2004, in prep.) to sulfur bearing species. SD, HDS, D 2 S, HDCS 
and D 2 CS, SD+, HDS+, D 2 S+, H 2 DS+, HD 2 S+ ,D 3 S+, HDCS+, D 2 CS+, H 2 DCS+, HD 2 CS+ 
and D 3 CS + have been added to the previous chemical network. We have assumed, as is usual 
in similar studies, that deuterated species behave as the main isotopic species and we have 
calculated the different branching ratios from statistical arguments except when dissociative 
recombination reactions were considered. In this latter case, we have assumed that removal 
of an hydrogen atom was more efficient that the release of D by a factor of 2. This has indeed 
been seen in various experimental studies of multiply deuterated molecular ions as discussed 
by Le Petit & Roueff (2003). The chemical reaction rate coefficients have been taken from our 
own chemical network (Pineau des Forets et al. 1993), updated with the recent compilation 
of Anicich (2003) and from the available chemical network of the Ohio State University 
(Herbst and coworkers). As is recalled in Millar (2003), deuterium enrichment is initiated 
by the fractionation reactions of H^, CH^ and C 2 H 2 f with D and HD. We have extended the 
reactions to the partially deuterated species until the fully deuterated ions D3 , CD3 and 
C 2 D 2 f are formed -as studied in the laboratory by Gerlich & Schlemmer (2002). Another 
important addition consists in considering the fractionation reactions with D 2 , which allows 



-13- 



more substantial deuterium enhancement. Recent laboratory studies have shown that the 
reactions are significantly less efficient as currently assumed (cf. Millar 2003) and that the 
reverse endothermic reaction may be more rapid than expected from pure thermodynamics. 
We have thus taken the reaction rate coefficients derived from the experiments (Gerlich 
et al. 2002) but determined the reverse reactions with the calculated endothermicity. We 
have also introduced possible fractionation reactions of H 2 CS + , H 3 CS + and H 3 S + with HD, 
which could occur but have not been considered until now as no experimental or theoretical 
information is available. The corresponding reactions and rate coefficients are displayed in 
Table 8 where the values of the endothermicities involved in the backward reactions have 
been guessed. As the Langevin collision rates have been assumed, the results obtained by 
including these reactions will give upper limits for the deuterium fractionation. 

Their influence is discussed below. We did not include the chemistry of methanol. 
The resulting chemical network consists in 210 species linked by more than 3000 chemical 
reactions. Details and discussion on the deuterated chemistry are discussed in Roueff et al. 
(2004, in prep.) and we focus the present discussion to the sulfur compound deuteration. 
Recent observations have shown that enhanced deuteration is linked to large depletions 
(see for example Bacmann et al. 2003). We consider then three different models with H 2 
densities of 10 4 , 10 5 and 10 6 and different elemental depletions as displayed in Table 9. These 
three models correspond approximately to the conditions found in TMC-1, Barnard 1 and 
L134N (see, Cernicharo & Guelin 1987; Bachiller et al. 1990, Pagani et al. 2004, in press). 
The cosmic ionization rate and the temperature are kept constant and equal respectively 
to 2xl0~ 17 s _1 and 10 K. The sulfur abundance is also taken constant and typical of the 
low metal case, i.e. S/H2 = 3.70xl0~ 7 . The Table shows that, as expected, deuteration 
increases from models 1 to 3. Fractionation reactions involving H2CO + , H3CO + , H3S + and 
H 3 CS + , are seen to be of considerable importance for the fractionation of H 2 S. However, 
their role is much less for formaldehyde and thioformaldehyde. This is readily explained 
by the chemistry as hydrogen sulfide is directly obtained from dissociative recombination of 
H 3 S + and involves a succession of ion molecule reactions. In the case of formaldehyde and 
thioformaldehyde, neutral reactions involving CH 3 and atomic oxygen and atomic sulfur are 
supposed to occur at a rate of 1.4xl0~ 10 . Then, deuterium enhancement of H 2 CO and H 2 CS 
reflects directly the deuterium enhancement of CH 3 . The reaction of CH 3 and O has only 
been studied in the laboratory down to 300 K and the reaction rate coefficient with sulfur 
is taken from the chemical network of the Ohio State University in the group of E. Herbst. 

As can be seen from the Table 9, this gas phase chemical model gives the correct orders 
of magnitude of the deuterium fractionation of ammonia, HCO + , N 2 H + for the high densities 
and large depletion factors when comparing to the observed values. In this model grains play 
a passive role, as they have trapped an important fraction of carbon and oxygen, allowing 
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the gas phase chemical processes to proceed. We also give the fractional abundances of 
formaldehyde, sulfur hydride and thioformaldehyde to check the relevance of our gas phase 
chemical network. We note that the fractional abundance of sulfur hydride is between 10~ 9 
and 4xl0 -10 for the three models. This compares reasonably to the values found in TMC-1 
and L134N as reported by Vastel et al. (2003). Grain surface chemistry is often invoked 
for the formation of H 2 S as the S + + H 2 reaction is highly endothermic. In the present 
chemical network, SH + is formed via S + Hg", which subsequently reacts with H 2 via radiative 
association to produce H 3 S + . H 2 S is then obtained from the dissociative recombination of 
H 3 S + . The fractional abundance derived for H 2 CS is below the observed ones by values of a 
few. However, we have checked that increasing the elemental abundance of sulfur results in 
increasing the steady state fractional abundances of H 2 S and H 2 CS. We do not try to obtain 
a perfect fit with the observed values as many physical and chemical parameters are not well 
constrained. The present stage of agreement is nevertheless encouraging. 



5. CONCLUSIONS 

We have carried out a comprehensive observational study in Barnard 1 of the deuterated 
species of thioformaldehyde, together with the observation of several lines from the main 
isotope and from the 34 S isotopic substituted species. In order to compare the deuterium 
enhancement we have also observed formaldehyde and its deuterated species. The main 
results are: 

1) The systematic velocity of the dense gas in Barnard 1 has been determined to be 
6.70±0.05 km s~ x from the observation of near 20 optically thin molecular lines with well 
known frequencies. 

2) The observed frequencies for HDCS, D 2 CS and H 2 C 34 S as well as the available lab- 
oratory data have been fit to obtain a set of molecular rotational constants allowing precise 
frequency determinations for these species in the millimeter domain. 

3) Rotational temperatures and column densities for H 2 CS, H 2 CO and their isotopes 
have been derived from rotational diagrams or through the use of a LVG code. The de- 
rived H 2 CS/HDCS and H 2 CS/D 2 CS abundance ratios are ~3 and 9 respectively, while for 
formaldehyde are ~7 and 13. 

4) Steady state gas phase chemical models including multiply deuterated species and 
with appropriate depletion factors, such as those found in different dense cores, allow to 
reproduce the deuterium fractionation approximately. In this context, surface processes 
have not to be included specifically , as they are poorly known, and their role is essentially 
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passive as the abundant oxygen, CO, H 2 ... components are frozen on the grains and are 
not available for gas phase chemical processes to proceed. Our results involve also some 
specific chemical assumptions which contribute to the high deuterium fractionation, namely 

(1) specific branching ratios in dissociative recombination of deuterated molecular ions with 
a 1:2 ratio for the ejection of H compared to D as observed in several experiments. And 

(2) fractionation reactions of H 2 CO + , H 3 CO + , H 3 S + , H 2 CS + with HD. These processes 
have not been studied in the laboratory nor theoretically and such studies would be very 
desirable. However, their influence on our results is not determinant except for the deuterium 
fractionation of H 2 S. 

We would like to thank Spanish Ministry of Science and Technology, MCyT, for funding 
support under grants ANAYA2000-1784, ANAYA2003-2785 and PNIE2001-4516. The au- 
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Fig. 1. — D 2 CS transition for Barnard 1 (dashed line) and for TMC-1 (solid line). The 
position of TMC-1 is that of cyanopolyyne peak RA(1950.0)= 04 h 38 m 38 s .O, dec(1950.0) = 
+25°35'45"; that of Barnard 1 is RA(2000.0)= 03 h 33 m 20 s .8, dec(2000.0)= +31°07'34". 



-20- 




Fig. 2. — Observed transitions for HDCS (dashed line) and D 2 CS with the frequencies 
originally adopted from the literature (a) and the corrected frequencies (b). 
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Fig. 3. — Observed transitions for H 2 CS. The velocity scale in the left panels (A) has 
been derived from the laboratory frequencies, while for the central panels (B) the predicted 
frequencies from a fit to the laboratory data were used (see text). In the plots of these two 
columns the horizontal lines represent the 3cr uncertainty for line frequencies (measured or 
predicted). Finally right panels (C) show the H2CS lines when a velocity of 6.7 km -1 is 
assumed for the main velocity component. 
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Fig. 4. — Observed transitions of H 2 C 34 S in the direction of Barnard 1. 
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Fig. 5. — Observed transitions from the different isotopic species of formaldehyde. Velocity 
scale has been derived for the main isotope from the calculated frequencies for this species, 
while for the rare species the measured space frequency, assuming a velocity of 6.7 km" 1 , 
have been adopted (see text). 
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Fig. 6. — Rotational diagrams for HDCS, D 2 CS and H 2 CS. The errors are given between 
brackets and concerns the last significative digit. Column density values are the averaged 
ones. 
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Table 1: Line parameters for the observed transitions of thioformaldehyde species 
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a Observed frequencies have been estimated from a Gaussian fit to the lines. 
Quoted frequencies correspond to the velocity component at 6.7 km s _1 
h Errors for observed frequencies are given between brackets 
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Table 2: Line parameters for the observed transitions of formaldehyde species 
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a Frequencies from Lovas catalog (see text) 
^Different line components 

c Complex line profiles with wings and selfabsorption 
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Table 3: Rotational constants for H 2 CS 



Parameter 


Laboratory 


Lab+Sky 


A 


291718.6 (133) 


291708.3 (169) 


B 


17699.6447 (146) 


17699.6529 (168) 


C 


16651.8711 (146) 


16651.8587 (165) 


Aj 


2.0638 (391) x 1(T 2 


1.8916 (178) x 10~ 2 


A JK 


5.20998 (739) x 10" 1 


5.22872 (573) x 10" 1 




22.8519 a 


22.8519 a 


8 3 


1.21460 (129) x 10~ 3 


1.21438 (163) x 10~ 3 




3.1922 (658) x 10" 1 


3.2438 (836) x 10" 1 




1.056 (492) x 10~ 5 


-8.41 (313) x 10" 6 



"Value adopted from Clouthier & Ramsay (1983) 
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Table 4: Rotational constants for D 2 CS 



Parameter 


Laboratory 


Lab+Sky 


A 


146400.549 (284) 


146397.685 (725) 


B 


14904.7104 (116) 


14904.70899 (556) 


C 


13495.4298 (116) 


13495.38776 (506) 


Aj 


1.2579 a x 10~ 2 


1.29708 (846) x 10~ 2 




2.8869 (166) x 10" 1 


2.7972 (178) x 10" 1 


A K 


5.2084 b 


6.588 (294) 


S 3 


1.40587 (269) x 10~ 3 


1.38673 (369) x 10~ 3 




2.14468 (275) x 10" 1 


2.2554 (245) x 10" 1 




1.198 (106) x 10~ 6 


7.25 (125) x 10- 7 




-3.6562 (919) x 10~ 4 


-3.075 (151) x 10~ 4 



"Adopted value from Cox et al. (1982) 
b Force field value (Turner et al. 1981) 
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Table 5: Rotational constants of HDCS 



Parameter 


Laboratory 


Lab+Sky 


A 


202716.66 (217) 


202716.54 (319) 


B 


16111.38339 (610) 


16111.38820 (887) 


C 


14890.97584 (603) 


14890.97943 (867) 


Aj 


1.587381 (547) x 10~ 2 


1.588513 (787) x 10~ 2 


A JK 


3.180421 (533) x 10" 1 


3.181521 (768) x 10" 1 




14.72 a 


14.72 a 


5 3 


1.38002 (484) x 10~ 3 


1.37782 (621) x 10~ 3 


h. 


2.8831 (281) x lO" 1 


2.8839 (414) x 10" 1 


Hj^ 


2.028 (217) x 10~ 6 


2.476 (312) x 10" 6 




-4.188 (109) x 10~ 5 


-4.193 (160) x 10~ 5 



a Intcrpolated value from H2CS and D2CS 
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Table 6: Rotational constants of H 2 C 34 S 



Parameter Laboratory Lab+Sky 

A 291708.3 a 291708.3 a 

B 17389.1451 (120) 17389.1477 (108) 

C 16376.6721 (120) 16376.6745 (108) 

Aj 1.8417 (486) x 10~ 2 1.8283 (399) x 10~ 2 

A JK 5.0236 (169) x 10" 1 5.0328 (138) x 10" 1 

A K 22.8519 a 22.8519 a 

8j 1.1662 (158) x 10~ 3 1.1673 (151) x 10~ 3 

6 k 3.2438 a x 10- 1 3.2438 a x 10" 1 



a Fixed to the value of H 2 CS 



-31 - 



Table 7: Relative abundances of deuterated molecules in several sources 



Molecule 


Barnard 1 


TMC-l a 


TMC-1 6 


L134N 


IRAS 16293-2422 


Orion (KL) 


DCO+/HCO+ 




0.012 1 




0.180 2 


0.009 3 


<0.01 4 


N 2 D + /N 2 H+ 


0.150 5 


0.0064 1 


0.080 2 


0.300 5 




<0.30 4 


NH 2 D/NH 3 


0.128 6 


<0.001 1 


0.020 2 


0.10 2 ' 7 


0.100 3 


0.062 8 


ND 2 H/NH 2 D 








0.050 7 






ND3/NH3 


0.0008 9 












HDCO/H 2 CO 


0.143 




0.059 1 


0.068 1 


0.140 3 


0.140 8 


D 2 CO/HDCO 


0.400 








<0.50 7 


0.021 8 


HDS/H 2 S 










0.100 3 


<0.02 n 


D 2 S/HDS 


0.108 12 








<0.264 12 




HDCS/H 2 CS 


0.300 


0.020 13 










D 2 CS/HDCS 


0.333 












CH3OD/CH3OH 
CH 2 DOH/CH 3 OH 
CHD 2 OH/CH 3 OH 
CD3OH/CH3OH 




0.026 1 




<0.032 1 


0.015 14 
0.306 14 
0.060 14 
~0.010 14 


0.01-0.06 15 



a Cyanopolyyne peak 
b Ammonia peak 

References. - 1. Turner (2001); 2. Tine ct al. (2000); 3. van Dishocck et al. (1995); 4. Turner & 
Zuckerman (1978); 5. Gerin et al. (2001); 6. Saito et al. (2000); 7. Roueff et al. (2000); 8. Turner (1990); 
9. Lis ct al. (2002); 10. Loinard et al. (2000); 11. Minh ct al. (1990); 12. Vastel et al. (2003); 13. Minowa 
et al. (1997); 14. Parise et al. (2004); 15. Mauersberger et al. (1988). 
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Table 8: Fractionation reactions introduced in the network 



Reaction 








P 




(cm 3 




(K) 


H 3 CO+ + HD ^ H 2 DCO+ + H 2 


2.0 


x 


10" 


-9 


500 


H 3 S+ + HD ^ H 2 DS+ + H 2 


2.0 


x 


10" 


-9 


500 


H 2 DS+ + HD ^ HD 2 S+ + H 2 


2.0 


x 


10" 


-9 


400 


HD 2 S+ + HD ^ D 3 S+ + H 2 


2.0 


x 


10" 


-9 


300 


H 3 S+ + D 2 ^ HD 2 S+ + H 2 


2.0 


x 


10" 


-9 


500 


H 2 DS+ + D 2 ^ D 3 S+ + H 2 


2.0 


x 


10" 


-9 


400 


H 2 DS+ D 2 -> HD 2 S+ + HD 


2.0 


x 


10" 


-9 


400 


HD 2 S+ + D 2 ^ D 3 S+ + HD 


2.0 


X 


10" 


-9 


300 


H 2 CS+ + HD ^ HDCS+ H 2 


2.0 


x 


10" 


-9 


500 


HDCS+ + HD ^ D 2 CS+ H 2 


2.0 


x 


10" 


-9 


400 


H 2 CS+ + D 2 ^ D 2 CS+ H 2 


2.0 


X 


10" 


-9 


500 


H 2 CS+ + D 2 ^ HDCS+ HD 


2.0 


x 


10" 


-9 


400 


HDCS+ + D 2 ^ D 2 CS + HD 


2.0 


x 


10" 


-9 


300 


H 3 CS+ + HD ^ H 2 DCS+ H 2 


2.0 


x 


10" 


-9 


500 


H 2 DCS+ + HD ^ HD 2 CS+ H 2 


2.0 


X 


10" 


-9 


400 


HD 2 CS+ + HD ^ D 3 CS+ H 2 


2.0 


X 


10" 


-9 


400 


H 3 CS+ + D 2 ^ HD 2 CS+ H 2 


2.0 


X 


10" 


-9 


500 


H 3 CS+ + D 2 ^ H 2 DCS+ H 2 


2.0 


X 


10" 


-9 


400 


H 2 DCS+ + D 2 ^ HD 2 CS+ HD 


2.0 


X 


10" 


-9 


500 


HD 2 CS+ + D 2 ^ D 3 CS + HD 


2.0 


X 


10" 


-9 


500 


H 3 S+ + D ^ H 2 DS+ + H 


2.0 


X 


10" 


-9 


500 


H 2 DS+ + D ^ HD 2 S+ + H 


2.0 


X 


10" 


-9 


400 


HD 2 S+ + D ^ D 3 S+ + H 


2.0 


X 


10" 


-9 


300 


H 2 CS+ + D ^ HDCS+ + H 


2.0 


X 


io- 


-9 


500 


HDCS+ + D ^ D 2 CS+ + H 


2.0 


X 


10" 


-9 


400 


H 3 CS+ + D^ H 2 DCS+ + H 


2.0 


X 


10" 


-9 


500 


H 2 DCS+ + D ^ HD 2 CS+ + H 


2.0 x 10" 


-9 


400 


HD 2 CS+ + D ^ D 3 CS+ + H 


2.0 


X 


10" 


-9 


300 



"The forward rate is kf, the reverse rate, k r , is given by k/ e ®I T 
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Table 9: Prediction of sulfur containing molecules fractional abundances and deuterium frac- 
tionation. 



model conditions 


Model 1 


Model 2 


Model 3 




H 2 (cm -3 ) 


10 4 


10 5 


10 6 




C/H 2 


7.5 x KT 5 


3.0 x 10~ 5 


1.5 x 10- 


•5 


0/H 2 


2 x KT 4 


8 x 10~ 5 


4 x 10~ 5 




N/H 2 


2 x 10- 5 


2 x 10~ 5 


2 x 10~ 5 




representative molecular cloud 


TMC-1 


Barnard 1 


L134N 




steady state results 


H (cm 3 ) 


2.24 


2.22 


2.25 




HD (cm- 3 ) 


2.76 x KT 1 


2.76 


26.8 




x(e") a 


6.03 x 10" 8 


3.52 x 10" 8 


2.60 x 10 


-8 


x(H 2 CO) 6 


1.2 x 10" 9 


4.7 x 10~ 10 


1.1 x 10- 


10 




1.3 x KT 9 


5.0 x 10~ 10 


1.1 x 10" 


•10 


x(H 2 S) 6 


4.5 x KT 10 


8.1 x 10~ 10 


8.7 x 10- 


•10 




6.1 x lO" 10 


1.4 x 10" 9 


1.7 x 10" 


9 


x(H 2 CS) 6 


4.1 x KT 11 


2.9 x 10- 11 


1.0 x 10- 


11 




4.8 x lO" 11 


2.9 x 10~ n 


1.3 x 10- 


•11 


D/H 


0.0053 


0.012 


0.021 




D 2 /HD 


0.021 


0.036 


0.059 




H 2 D+/H 3 + 


0.029 


0.10 


0.16 




D 2 H+/H 2 D+ 


0.054 


0.10 


0.16 




D 3 + /D 2 H+ 


0.041 


0.086 


0.15 




HDCO/H 2 CO fe 


0.091 


0.12 


0.14 






0.048 


0.07 


0.089 




D 2 CO/HDC0 6 


0.34 


0.57 


0.80 






0.60 


0.95 


1.20 




HDS/H 2 S 6 


0.26 


0.37 


0.43 






0.017 


0.032 


0.050 




D 2 S/HDS 6 


0.53 


1.0 


1.32 






0.0049 


0.007 


0.0086 




HDCS/H 2 CS 6 


0.16 


0.18 


0.18 






0.041 


0.065 


0.086 




D 2 CS/HDCS 6 


0.48 


0.85 


1.17 






0.37 


0.67 


0.92 




DCO+/HCO+ 


0.028 


0.043 


0.066 




NH 2 D/NH 3 


0.040 


0.077 


0.20 




ND 2 H/NH 2 D 


0.027 


0.041 


0.051 




ND 3 /ND 2 H 


0.035 


0.05 


0.056 




N 2 D+/N 2 H+ 


0.032 


0.047 


0.068 





